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High-Energy Gamma-ray Emission Properties of an Anomalous X-ray Pulsar, 

4U 0142+61 
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Sabanci University, FENS, Orhanli- Tuzla, 34956, Istanbul, Turkey 

Anomalous X-ray Pulsars (AXPs) are bright X-ray sources. Few AXPs emit highly pulsed emission 
in hard X-rays. Using data collected with the Large Area Telescope on board Fermi Gamma-ray 
Space Telescope, we explored high-energy gamma-ray emission from the brightest AXP, 4U 0142-1-61. 
We do not detect any significant emission from the source. Here, we present the upper limits to the 
persistent and pulsed emission of 4U 0142-1-61 in the high-energy gamma-ray domain. 
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I. INTRODUCTION 

Anomalous X-ray Pulsars (AXPs) are bright X-ray 
emitters; their luminosities (below 10 keV) range from 
lO'^'^ to 10'^^ erg s~^ that exceed their spin-down lumi- 
nosity. Their spin periods are clustered between 2—12 
s and spin-down rates are relatively large, i.e., between 



10"^'^ — 10"^"^ s (see [l| for a review). Their obser- 
vationalproperties are explained either with a fallback 
disk [5, 6\ or by the decay of their strong magnetic field 
(magnetar model) Q. 

Until the discovery of hard X-ray emission from 
AXPs IE 1841-045 [IL 4U 0142-^61 Q and 
IRXS J170849.0-400910 Q they were known to be 
emitting soft X-rays only (< 10 keV). The origin of 
hard X-ray emission is still not well understood. Emis- 
sion from corona [8j], breakdown of fast- modes via 
quantum electrodynamics effects d, [l0[ and resonant 
Compton upscattering of soft photons [TTI| are pro- 
posed to explain the hard-X ray emission. 

Construction of spectral behaviors of these sources 
in a broad energy band is important for understand- 
ing the origin of the hard X-ray emission. Here we 
present the upper limits to the persistent and pulsed 
high-energy gamma-ray emission from 4U 0142-1-61 at 
the GeV range. We estimated an upper limit to the 
spectral break energy by extending the vF^, spectrum 
of the source to GeV range. 



II. OBSERVATIONS AND DATA ANALYSIS 

The LAT observations between 2008 August 4 to 
2010 April 29 with an exposure time of ~31.7 Ms 
were used to investigate the persistent emission from 
4U 0142+61. We obtained and analyzed the data from 
15° radius around the source and also, a 2° radius 
region was selected and analyzed in order to avoid 
contamination from nearby bright sources. Spectral 
fits and flux calculations were done with the Python 
version of gtlike, pyLikelihood, for the 0.2—1.0 and 
1.0-10.0 GeV energy bands. See [2l| for the details 
of LAT data calibration. 

To perform timing analysis, all event photon arrival 
times were extracted from a 2° region and converted 



by gtbary to the arrival times at the solar system 
barycenter. We used 2—10 keV RXTE observations 
and found the spin ephemeris of the source between 
2008 August 4 and 2010 April 30 which also covers 
the extracted LAT observations. 

To search for pulsed high-energy gamma-ray emis- 
sion from 4U 0142-1-61, first we obtained the spin 
ephemeris of the source using contemporaneous 
RXTE/PCA observations in the 2-10 keV range with 
a total exposure of 196 ks. A Fourier based epoch fold- 
ing technique was applied to the data to obtain the 
spin ephemeris. We generated the pulse profiles of 
the source using three consecutive PC A observations 
around the epoch (MJD 54713.5) which are grouped 
such as they are separated at least 0.2 days from each 
other. We determined the phase shift of each pointing 
with respect to the template by cross correlating the 
pulse profiles of each group of pointings with the tem- 
plate profile and fitted the phase shifts with a polyno- 
mial. In Table 1 we present the best fit spin ephemeris 
parameters of 4U 0142-1-61. We used the precise PCA 
spin ephemeris that we obtained to search for pulsed 
high-energy gamma-ray emission from 4U 0142-1-61. 
We generated the LAT pulse profiles in the 0.2—1.0 
GeV and 1.0—10.0 GeV energy ranges and found that 
both LAT profiles are consistent with random fiuctu- 
ations with respect to its mean. 



TABLE I: Spin ephemeris of 4U 0142-f61 as determined 
using RXTE/PCA observations. 



Parameter 



Value 



Range (MJD) 
Epoch (MJD) 
V (Hz) 

u (10-'^ Hz s-^) 
j> (10-^^ Hz s"^) 



54682.6 - 55315.1 
54713.5 
0.1150900026(9) 
-2.745(8) 
3.6(3) 



III. RESULTS 

After processing the data as explained in §2, we 
fitted a power-law to the data obtained from 15° 
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FIG. 1: A wide band spectrum of 4U 0142+61: INTEGRAL /ISGRl (20-300 keV) in black (stars), INTEGRAL/SPI 
(20-1000 keV) in red (open squares) and CGi? 0/COMPTEL (0.75-30 MeV) 2a upper limits in black (data obtained 
from [13]). Blue diamonds are the Fermi/LAT upper limits in the 0.2—1.0 GeV and 1.0—10.0 GeV obtained using 2° 
extraction region. Red triangles are upper limits for the 15° extraction region. Dashed line is the best fit power-law 
model to the ISGRI data points [l^]. Solid line shows the power- law upper limit trend of the 2° Fermi /L AT region. 
Figure is from [2l[ 



radius region with an index of 2.5. The fit results a 
test statistics (TS) value of ~0.23 which implies a 
detection significance less than la. We calculated 3a 
flux upper Hmits as 2.32 x 10~^ MeV cm~^ s~^ in the 
0.2-1.0 GeV band and 1.28 x IQ-^ MeV cm-^ 
in the 1.0—10.0 GeV band. For the 2° radius region 
power-law fit with an index of 3 resulted in a TS 
value of '-.^3 which implies a detection significance less 
than 2a. The 3a flux upper limits are 5.72 x 10~^ 
MeV cm-2 s'^ and 1.29 x 10"^ MeV cm'^ s^^ for 
0.2-1.0 GeV and 1.0-10.0 GeV, respectively. 

We searched for pulsed high-energy gamma-ray 
emission from 4U 0142-1-61 by folding the Fermi/LAT 
data with the precise spin ephemeris obtained using 
contemporaneous RXTE/PCA observations. Our 
search for emission yields random fluctuations with 
respect to its mean (See Figure [2]). The 3a upper 
limits to the RMS pulsed amplitude are 1.5% and 
2.3% in the 0.2-1.0 GeV and 1.0-10.0 GeV band, 
respectively. A search in the lower energy part of 
the LAT passband (30-200 MeV) also yields in no 
evidence of pulsed emission; the 3a RMS pulsed 
amplitude upper limit is 1.6%. 



IV. DISCUSSION 

We used the LAT persistent emission upper limits 
to extend the vFi, spectrum of 4U 0142-1-61 presented 
in [l3] to the 10 GeV. We could place an upper limit to 
the very high-energy gamma-ray emission trend (i.e., 
the LAT results obtained using the 2° region) which 
is a power law with an index of —0.76 (solid line in 
Figure [I]). INTEGRAL/lSGRl data is fitted with a 
power law index of 0.93 ± 0.06 by ^1?| as shown with 
dashed line in Figure[TJ The intersection of two curves 
suggests an upper limit of '^l.l MeV to the spectral 
break energy. 

The models presented in Q and [l^l estimate a 
spectral break energy of about 1 MeV which is con- 
sistent with the break energy that we estimated. If 
there is a significant excess emission in optical band 
as in the 4U 0142-1-61 [l^, quantum electrodynamics 
model [l,[l3| predicts an increase in the 10 — 200 MeV 
range. Excess optical emission can be due to the star 
itself [11 [Hi or the disk around the star [H, [i3|- If 
the source of the excess optical emission is the disk, 
quantum electrodynamics model is consistent with our 
results since the star itself would not be excessive in 
optical band. Outer gap model of AXPs within the 
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FIG. 2: Pulse Profiles of 4U 0142+61 by folding the Fermi/LAT data with the precise spin ephemeris obtained with 
contemporaneous RXTE/PCA observations. Left: in the 30 - 200 MeV band, middle: in the 200 MeV - 1 GeV band, 
Right: in the 1-10 GeV band. 



magnetar interpretation ^ also predicts high-energy (MTKD-CT-2006-042722). 
gamma-ray emission from the source. Recently, [19| 
and [1^ suggested that non-detection is inconsistent 
with the outer gap model and fallback interpretation 
cannot be eliminated. 
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